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. Abstract. We review Migdal's Theory of Finite Fermi Systems and its application to the 

\ structure of nuclei. The theory is an extension of Landau's Theory of Interacting Fermi 

^> ■ Systems. In the first part the basic formulas are derived within the many body Green functions 

approach. The theory is applied to isovector electric giant resonances in medium and heavy 
' mass nuclei. The parameterizations of the renormalized effective ph-interaction and the 

effective operators are discussed. It is shown that the number of free parameters are restricted 
" ' \ due to conservation laws. We also present an extension of Migdal's theory, where the low- 

■ lying phonons are considered in a consistent manner. The extended theory is again applied to 
l/^ , the same isovector electric giant resonances and to the analysis of (oc,Oc') reaction data. We 

' point out that the extended theory is the appropriate frame for self consistent nuclear structure 

' calculations starting from effective Lagrangians and Hamiltonians. 

\o '■ 
o ■ 

■ 

^ i 1. Introduction 

^ ■ The Landau-Migdal theory of Finite Fermi Systems^^ is the best- founded phenomenolog- 
' ical theory for nuclear structure calculations. It is based on Landau's Theory of Interacting 
is> ' Fermi Systems ||3l. The basic formulas can be derived within the many-body Green func- 

■ tions formalism in an exact way using Landau's renormalization procedure. This procedure 
. give rise to renormalized residual interactions and effective single particle operators. The 

basic input which is required to solve the corresponding equations are the quasi particle ener- 
gies, the renormalized particle-hole(ph) interaction and effective operators. If one considers 
pairing-correlations also a renormalized particle-particle interactions enters. So far in most of 
the applications the quasi particle energies have been taken from experiment and/or single- 
particle models and the renormalized residual interactions have been parametrized and the 
corresponding parameters are fitted to data. These Landau-Migdal parameters turned out to 
be universal for all nuclei, except the lightest ones. More over some of these parameters 
are directly connected with gross properties of nuclei such as the compression modulus and 
the symmetry energy. For that reason any effective interaction, Hamiltonian or Lagrangian 
which is used for nuclear mean field calculations should be tested against specific Landau- 
Migdal parameters. This possibility was first used in 1975 by Backman, Jackson and Speth 
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Mi who showed that several parameterizations of the Skyrme forces used at that time led to 
an instability against spin collapse and had to be modified. In Ref. this method was used 
to design new versions of the Skyrme interaction which could be used in self consistent nu- 
clear structure calculations. In Ref. ||6| the Landau-Migdal parameters where calculated from 
G-matrices based on meson exchange models. The authors showed in addition that the Lan- 
dau expansion of the ph-interaction is well justified for all parameters except the spin-isospin 
parameter G' which is strongly momentum dependent due to the long range of the one-pion 
exchange. For that reason one has to include this contribution explicitly if one calculates mag- 
netic properties iQ. A detailed discussion of the Landau Migdal theory and self consistent 
approaches and the application of effective forces in the FFS can be found in the review by 
Khodel and Saperstein lISj. Fully self consistent calculations in the Landau-Migdal approach 
have been also performed by Fayans et al. [91 who based their mean field and the residual 
interactions on local energy -density functionals. Chiral Lagrangians lITOll seems to be espe- 
cially appropriate for self consistent nuclear structure calculations within the Landau-Migdal 
approach because they have a similar structure as the generalized Landau-Migdal interaction 
(which includes the one pion exchange contribution [6]). They have strong short range com- 
ponents and the one pion exchange is included explicitly. It has been shown in |6| that this 
structure of the bare interaction does not change if one performs a partial summation to obtain 
the Brueckner G-matrix. The inclusion of the one pion exchange in chiral Lagrangians is the 
basic difference to Skyrme type Hamiltonians which are commonly used in self consistent 
nuclear structure calculations. This ansatz is a parametrization of a G-matrix and includes 
zero-range and derivatives of zero range components only. The zero range three body part can 
be replaced by a density dependent two-body interaction which simulates the density depen- 
dence of the interaction in nuclei. In the case of chiral Lagrangian the density dependence is 
introduced by the G-matrix and more over it has been shown in Ref.[ 10 1 that the chiral expan- 
sion also give rise to a strong three body force if applied to the nuclear many-body problem. 
It is important to mention that the single particle spectrum calculated in a mean field approach 
can not be directly compared with the experimental spectrum and can therefore not be iden- 
tified with Landau's quasi particles because the coupling to phonons modifies the theoretical 
spectrum. For that reason, as we will show, the mean-field spectra are the appropriate input 
for the Extended Theory of Finite Fermi Systems [1 IJ, which includes in a consistent way the 
phonon contributions. 

2. Landau's Basic Theory 

Landau's theory of interacting Fermi systems Q deals with weakly excited states in infinite 
systems of interacting fermions, such as ^He or nuclear matter. Without the interaction the 
system would be simply a collection of independent fermions, each characterized by its 
momentum p and by its spin. The particles obey the Fermi-Dirac statistics and therefore 
occupy all quantum states up to the Fermi momentum pf \ all the states with larger momenta 
are unoccupied. The occupied states form a sphere of radius pp, which is called the Fermi 
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sphere. The Fermi momentum is proportional to the number of particles 



?,n^h^N 



1/3 



Excited states are formed by moving particles from the interior to the exterior of the Fermi 
sphere. Therefore excited states differ from the ground state by a number of particle-hole 
pairs, the particles with p > Pf and the holes with p < pf. These particles and holes play the 
role of elementary excitations in the ideal Fermi gas. They appear and disappear only in pairs 
and they have momenta close to pp for weakly excited states. 

The basic assumption by Landau is that interacting systems may be obtained from non- 
interacting systems by an adiabatic switching on of the interaction. In particular, he assumes 
a one-to-one correspondence between the single-particle states of the non-interacting system 
and the so-called quasi-particle states in the interacting system. These quasi-particles behave 
in the correlated system like real particles in the non-interacting system: They obey the Fermi- 
Dirac statistics and therefore occupy the corresponding quasi-particle states up to the Fermi 
momentum which can be calculated from the same formula (Eq.©). As in the non-interacting 
system, the excitations in the correlated system are (quasi-)particles with momenta greater 
than Pf and (quasi-)holes with momenta less than pf that appear and disappear, again in 
pairs. Quasi-particles are only defined near the Fermi surface, where their decay width is much 
smaller then their excitation energy. In order to define such quasi-particles Landau made the 
following assumption concerning the interaction of the quasi-particles: the interaction follows 
from a self-consistent field that acts on one specific quasi-particle and which is produced by 
all the other quasi-particles. For that reason the energy of the system is no longer the sum 
of the separate quasi-particle energies, but is a functional of the occupation function n{p) 
of the quasi-particle state p. If one excites the system one basically changes the occupation 
functions n{p) by an amount 5n(p). The corresponding change of the energy is 

5£ = j dp Z^{p)hn{p) + JdpJ dp'F"'{p,p')dn{p)dn{p') (2) 

where e^{p) are the equilibrium energies of the quasi-particles. We rewrite this expression as 

8E = j dp £°(p)+ j dp'F'^{p,p')8n{p') bn{p) = j dpe{p)8n{p) (3) 

The natural definition of the energy of an individual quasi-particle is therefore the first 
variational derivative of the total energy with respect to the occupation function and that of 
the interaction between the quasi-particles is the second functional derivative 

on{p) on{p)on{p') 
In ^He the interaction F^ has a spin-independent and a spin-dependent part. 

F"'{p,p')=F{p,p')+G{p,p')o-a' (5) 

The energies of the elementary excitations are measured from the Fermi sphere. For quasi- 
particles one obtains 
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and for quasi-holes 



The velocity of the excitations at the Fermi surface is given by 

PF 



(8) 



As we shall see later the effective mass m* is directly connected with F'^. 

If one considers excitations near the Fermi surface one may restrict the magnitude of the 
momenta p and p' in the interaction F^{p,p') to the Fermi momentum. The only variable 
then is the angle between the two momenta. This suggests an expansion in Legendre 
polynomials 

f{p,p')=Y^fiPi{cos&) (9) 

/ 

The constants // are the famous Landau parameters. One introduces dimensionless parameters 
by defining 



Fi=Cofi (10) 



where 



is the inverse density of states at the Fermi surface. Likewise, one may also expand the other 
component of the interaction in eq.® and define the parameter gi in the same way. The 
parameters fo and /i are related to the compression modulus K by 

r-,2 dHE/A) _ h^pl I .r. .... 
^-P^^p^-^^^^ + ^f^^ (12) 



and the effective mass m 

mVm =l + ^/i. (13) 

The interacting quasi-particles and quasi-holes in ^He give rise to the famous zero sound and 
spin waves. The corresponding linearized time dependent equation for dn{p) is equivalent to 
the random phase approximation (RPA) equation; it is the basis of the Landau- Migdal theory. 



3. Migdal's Theory of Finite Fermi Systems 

Migdal has applied Landau's theory of interacting Fermi system to atomic nuclei Q- Here 
one deals with two kinds of fermions -protons and neutrons- and, more importantly, with a 
relatively small number of particles. In the past 40 years Migdal's theory has been successfully 
applied to numerous nuclear structure problems. Whereas in infinite systems the quasi- 
particles differ from the bare particles only by their effective mass, the situation in finite 
nuclei is completely different. Here the quasi-particles are the single particle excitations of 
the nuclei, which can be obtained from phenomenological shell models or, in the case of 
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self consistent approaches, by the corresponding mean-field solutions of a given effective 
Lagrangian, effective Hamiltonian and energy density functionals, respectively. The quasi- 
particle interaction is defined in the same way as before. The interaction, however, depends 
not only on the spin but also on the isospin of the quasi-particles and is density dependent 
because the interaction inside a nucleus is very different from the interaction in free space. 
As in infinite systems, one expands the interaction at the Fermi surface in terms of Legendre 
polynomials and the parameters of the expansion -the famous Landau-Migdal parameters- 
are determined from experiment. These parameters turn out to be universal for all nuclei but 
the very light ones. 

As pointed out in the introduction,the Theory of Finite Fermi Systems (TFFS) is a semi- 
phenomenological microscopic theory. It is microscopic because all its equations are derived 
rigorously from first principles, however it contains also phenomenological aspects, such as 
the quasi-particle energies and wave functions and the quasi-particle-quasi-hole interaction. 
All these quantities are defined microscopically and could, in principle, be calculated starting 
from the bare nucleon-nucleon interaction. Such calculations are very involved and in actual 
calculations one has to make severe approximations, so that we cannot expect to obtain 
quantitative agreement between, for example, the calculated and phenomenological Landau- 
Migdal parameters. Nevertheless the calculated Landau-Migdal parameters are in surprisingly 
good agreement with the phenomenological ones 0. In the TFFS, as in Landau's theory, one 
derives a renormalized equation for the response functions and solves this equation in a given 
approximation. For that reason the TFFS can be applied to such problems as are connected 
with the response function. These include, for example, excitation energies and transition 
probabilities in even nuclei and moments and transitions in odd-mass nuclei. The input data, 
single particle energies and wave functions are supposed to be given by a shell model in the 
phenomenological approach, or are calculated in mean field approximation in self consistent 
calculations. Similarly the residual particle-hole and particle -particle interactions are either 
parametrized and adjusted to the experiment or determined from the corresponding effective 
Lagrangian and effective Hamiltonian, respectively. In this respect is the TFFS and its 
extension (which includes phonons explicitly) a frame for any nuclear structure calculations, 
for self consistent and non self consistent ones. This will be discussed in detail in the following 
chapters. 

3.1. Basic properties of many body Green functions 

In this section we will give a short introduction to the theory of many -body Green functions. 
(More details are given in |2^). They are the basis of the TFFS. The starting point of the theory 
is a Hamiltonian for Fermions 

// = ^£ya^av + — ^ l^V2V4ViV3 ^V4'^V2^Vi'2v3 (14) 
V ^Vi...V4 

Here are single particle energies, while ay and a^ are creation and annihilation operators 
of nucleons in the state cpv^ which obey anticommutation rules, and w is the residual two-body 
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interaction. We now consider an A-particle system 

H\An)=E\\An) (15) 

where \An) and E\ are assumed to be the exact solutions of the many-body Schrodinger 
equation with a given many-body Hamiltonian. The one-particle and two-particle Green 
functions are defined as 

gv,v2 {h;t2) = i-i) (A0|r{av,(?i)a+(^2)} |A0) (16) 

^viV3V2V4 {ht3;t2t4) = i-i)^) {A0\ T{ay^{ti)ay^{t3)c4^{t4)a+^{t2)} \A0) (17) 
with the time-dependent creation and annihilation operators 

av(f) = e'^' ay e''"' a+{t) = e'"' a+ e''^' (18) 

The symbol T denotes the time-ordering operator of Wick, which means that the operators 
should be taken in time-ordered form with the latest time to the left and the earliest to the 
right. In addition, the product has to be multiplied by a factor ±1 depending on whether an 
even or odd number of permutations is needed to come from the given to the time-ordered 
form. The Fourier transform of the one-particle Green function has the form 

1 r -\-oo 

gv,v2 ih;t2) = / rfa)^v,v2(o)) (19) 



271, 

Correspondingly, the two-particle Green function depends on three frequencies. For further 
investigation it is convenient to introduce Green functions with a source field q 

n o^ _ i-i){A0\T{Ua,,{h)a+{t2)}\A0) 
' ^ " {AO\T{U}\AO) ^ ^ 



with 



jj ^ J d5d6 q{5,6)a^^{ti)av^{tf,) ^2^) 



and the short hand notation: 

£ /'jfi= / dl. (22) 

In order to determine the Green functions one needs the equations of motion. They can 
be obtained by differentiation of the corresponding Green function with respect to time. For 
the one-particle Green function with a source field we obtain 

J J3{S(1,3) + q{\,3)}gq{3,2)+'-j J3J4J5w(14,35) ^^(35,24) =5(1,2) (23) 
with the abbreviations 

5v:V25(?i-?2)|^J^-eO,| =5(1,2). (24) 

Here we have used eq. (fT4ll in order to eliminate the Hamiltonian H in eq. (l23t . The equation 
for the two-body Green function is connected with the three-body Green function and so 
forth. Thereby one obtains a hierarchy of equations. In order to calculate the one-body Green 
functions one needs the two-body Green function, etc. This is not the way that we proceed 
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in the following; rather, we replace the bare two-body interaction w by an effective one-body 
potential E {mass operator or self-energy) which, in principle, is given by the bare interaction 
and the two-particle Green function. The corresponding equation for the one-particle Green 
function is called the Dyson equation. 

^|c/3 (5(1,3) +^(1,3) +E,(1,3)}^,(3,2) =5(1,2) (25) 

So far we have only replaced the unknown two-particle Green function by the also unknown 
self-energy E and little seems to be gained. Before we discuss the solution of this problem we 
continue with the development of the formal theory. We define the response function L by the 
functional derivative of the one-particle Green function 

5^,(1,2) 



L(13,24) = ^(13,24) -^(1,2)^(3,4) 



(26) 



6^(4,3) 

The functional derivative of the Dyson equation gives an integral equation for the response 
function 

(27) 



(28) 



L(13,24) = -^(1,4)^(3,2) -/y d5d6dld% ^(l,5)if(57,68)L(83, 74)^(6,2) 

where we introduced the effective two-body interaction K via 
5E(1,2) 



5^(3,4) 



/i^(13,24) 




K 



Figure 1. Graphical representation of the Bethe-Salpeter equation for the response function 
in the ph channel. 

A graphical representation of eq. (l?7t is given in FigCl One notices that the kernel K of 
the integral equation is irreducible with respect to the particle-hole propagator. We further 
introduce the vertex function x by 

L(13,24) = j d5d6 g{l,5)x{53M) g{6,2) 

The corresponding equation follows from eq. (l27b 

x(13,24) = -6(1,4)5(3,2) -^y" J5J6J7J8i^( 17, 28)^(8, 5)g(6,7)x(53, 64) 

The linear response function allows one to determine the following quantities : 

(I) From eq. (l26b we obtain in linear approximation the change of the one-particle Green 
function : 

5g,(l,2) = £l(13,24)5^(4,3) (31) 



(29) 



(30) 



3,4 



The Fourier transformation of the time variables gives dgq{QO,Q.) and, after integration 
over (0, one obtains Spq{Q.), where Q. is the energy transferred by the external field 
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dq{Q.) and dpq{Q.) is the change of the density due to this external field in a given 
single particle basis cpv '■ 

dp,,,,iQ.) = £ L(a)v 1V3V2V4 5<?V4V3(^^) (32) 

V3,V4 

(II) With the time order ^3,^4 > tij2, we insert into the two particle Green function a 
complete system of eigenfunctions of the A-particle system between the two particle- 
hole pairs 

00 

g(13,24) = £ /"(3%"0(12) (33) 

n=0 

with 

g0«(34) = -i{A0\T{a{3)a+{4)}\An) (34) 

and 

^"0(12) = -i{An\T{a{l)a+{2)}\A0) (35) 

As the response function is defined as 

L(13,24) =g(13,24) -g(34)g(12) (36) 

its spectral representation is also given by eq. (l33t but the sum begins atn= 1 . With the 
appropriate time order one obtains from the generalized Green function in eq. (l35l) the 
transition amplitudes between the ground state and the excited states of the A-particle 
system. 

(Ill) With the time order ±00 and t4 — > =1=0° in the two particle Green function one can 
move these operators to left and right side of the time ordering and may introduce two 
complete sets of eigenstates of the (A-i-1)- system or (A-l)-system, respectively. From 
this one obtains another generalized Green function 

g^fi(l2) = -i{A±l,a\T{a{l)a+{2)}\A±l,^). (37) 

where a and (3 are two given states in the odd mass system. As we know only the 
response function, one has to subtract the one-particle Green function of the A-particle 
system and obtains finally the change of the density between two given states in the 
(A ± l)-system and the ground state of the A-particle system. With this quantity one 
can calculate moments and transitions in odd-mass nuclei. Details are given in L2J. 

3.2. One-particle Green function and the nuclear shell model 

In order to solve the equation for the response function and vertex function, respectively, we 
have to know the one-particle Green functions ^(1,2) and the effective two-body interaction 
K. As Landau has shown, one does not need the full information included in g but only the 
so-called pole part of the one-body Green functions near the Fermi surface. The rest gives 
rise to a renormalized interaction and effective one-body operators. We first discuss the pole 
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part of the one-body Green function. For that we consider the Dyson equation in coordinate 
space Q 

(e- 1^) sikiXi.^) - J S(^i,^3,e) ^fe,^2,e) = 5(^i -^2) (38) 

where ^ = {r,s) represents space and spin coordinates. In an arbitrary single-particle basis 
(Pv(^) the Dyson equation has the form 

.2 



V3 



2m 



gVj,V2 — ^ViV2 



V1V3 



We now chose a special basis cpv(^) that diagonalizes the expression in the brackets 
.2 



2m 



(£)5v[V3 



(39) 



(40) 



V1V3 



It is obvious that such a basis has to depend on the energy e. The one-body Green function 
becomes diagonal in this basis and can be written as 

SviV2 



The functions cpv that diagonalize g satisfy the equation 

1 

|^-fE(^,p,£) (pv(^,e)=£v(e)cpv(^,e) 



(41) 



(42) 



Here the self-energy E is given in a "mixed" representation, where we have replaced the non- 
locality by a momentum dependence 



E(^,p,£) = J J^iE(^,^i,£) 



(43) 



Due to Landau's renormalization procedure one needs the one-particle Green functions only 
near the poles £v, with 

,2 



2m 



(44) 



vv 



The single-particle energies £v so defined are the quasi-particle energies (in the sense of 
Landau) for finite systems. We expand eq. (l42b around the quasi-particle energy £v up to 
first order in £ and obtain for the one-(quasi-)particle Green function 



^viV2(e) 



JV1V2 



1 



£-£vi+/Yvi fi 



(45) 



The residue of g at the pole is called the single particle strength 

1 

Zv = 



1 



dE 



(46) 



As we expanded the energy up to first order around the pole, one needs the wave functions 
only at the pole energy, (Pv(^) = 9v(^,ev)- These wave functions satisfy the equation 

|^-fE(^,p,£v) (pv(^)=evcpv(^) (47) 
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We introduce the quantity pj^, which is defined by 

E^+Z{r,pF{r),£F)=eF (48) 

where ep is the Fermi energy. In infinite Fermi systems pf can be identified with the 
momentum at the Fermi surface. We now expand the self-energy E about the Fermi 
momentum pp and the Fermi energy £f 

This formula can be rewritten as 

/ + [/(r) l(pv = evcpv (50) 



1 ,2 



2m* 

with the abbreviations 



m 



^ 5e 



(51) 



F 



and 



C/W=A'-4^ = -(SV+A^fl--) (52) 

> ^ 2m* m* ^ m*J 

If we also take the spin degree of freedom into account then, up to order in the self-energy 
E, one obtains the following single-particle Hamiltonian 

H = -^p^ + Uir)+ a(r)a ■ / + p(r)Z(Z + 1) (53) 
2m* 

Here 1 is the angular momentum operator and Z its quantum numbers. The form of this 
Hamiltonian follows from the expansion of the self-energy E about the Fermi energy ef 
and the Fermi momentum pf, and symmetry arguments. This Hamiltonian agrees with 
phenomenological shell model Hamiltonians. It is useful to point out that the spin-orbit term 
follows here from the expansion of the non-locality of E. 

3.3. Landau's renormalization procedure 

We now return to the one-particle Green function. With eqs. (I45landl46b we can separate the 
one-particle Green function into a singular part and a remainder 

^v.v. («) = ^ ^^r'^Vp T + ^v^v. (CO) (54) 



with 



Zvi = |(AO|avi|A±l vi )p 

£v, =<+l-£^ /or £vi>A^ (55) 
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With this ansatz one writes the product of two one-particle Green functions -a quantity that 
enters in all the integral equations we have derived before- as a singular part S and the 
remainder R, using the ansatz in eq.(l54b 

^ViV3 {^'^\^ •?V2V4 {^~\^ = 5'viV2V3V4 {^\^) +^ViV2V3V4 (56) 



with 



5'viV2V3V4 (ra;^^) = 27i:/Zv32V4SviV35v2V4 S ( (Q - + gv4 \ (57) 

\ ^ / ^V2 ~ ^Vi ~ 

Using Landau's renormalization procedure one obtains an equation for the renormalized 
vertex function x 

XviV3,V2V4 (to, ^) = 'CviV3,V2V4 ^) + (58) 



52 '^V^V5,V2V6 ( t 



The renormalized vertex is defined as: 

'CviV3,V2V4 (^5 ^) = \^ZviZv2 '^ViV3,V2V4 (^7 ^) (59) 

Here only the singular part of the product of the two Green functions appears explicitly, 
whereas the remainder gives rise to a renormalized effective two-body interaction F^^ and 
a renormalized inhomogeneous term x'^. In similar fashion one can renormalize the response 
function L and obtains 

L(13,24) =x''(13,57)L(57,68)x'^(68,24)+x'^(13,57)i?(57,24) (60) 

where L is the quasi-particle response function 

Z(13,24) =5(l,2)x(1324) (61) 

3.4. The renormalized equations of the TFFS 

As we have seen in section 2.1, the response function includes the transition amplitudes 
between the ground state and the excited states of an A-particle system. Using the arguments 
in section 2. 1 and the projection procedure described in pj, one obtains from the renormalized 
response function in eq. (l60b and the renormalized vertex function eq. (l58t an equation for the 
quasi-particle quasi-hole transition matrix elements 

(Evi — ev2 — ^) XviV2 = ("Vl ~^'^2) 52 '^viv4V2V3 XV3V4 (62) 

V3V4 

which are connected with the full ph-transition matrix elements 

%;;'0^=<Am|a+av2|A0> (63) 



by the relation 



XviV2 ^ '^V3ViV4V2 XV3V4 i^^) 
V3V4 
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The transition matrix element of a one-body operator Q is given by 

<Am|e|AO>= (65) 
V1V2 

where Q^^^ is the renormalized one-particle operator 

— H t;yjy^y2V42v3V4 (66) 
V3V4 

In section 2. 1 we also showed that the response function allows one to calculate the change of 
the density in an external field. 

5pViV2(^^) = LviV2V3V4(^^) S<?V3V4(^^) (67) 

V3,V4 

In the following we are interested in the transition strength as a function of the excitation 
energy. The expression for the strength function is given by 

S(t2,A) = --/m £ eviV2(^^) 5pv2Vi(^^ + iA) (68) 
^ V1V2 

where A is a smearing parameter that simulates the finite experimental resolution and the more 
complex configurations not included in the present approach. From eqs. (l58t to (l6TT) and (l66b 
we obtain, with dq = Q and 5p = p, 

S{a,A) = --ImY, Ql%{^)Pv,v,{n + iA) (69) 
^ V1V2 

where p is the change of the quasi-particle density 

PV,V2(^^) = /'\"'\ Ql%M - I <3,,,,Pv,V3(^^)(70) 

tvi— tv2— *^Vi ^ t-V2 — i"i V3V4 

The second term on the right side of eq. (l60b can be neglected because we are interested in the 
singular solutions (resonances) of the equation. The calculations presented here have been 
performed in r-space because this allows treatment of the continuum in the most efficient 
way. The corresponding equations have the form ifTTlfm 



p(r,a) =- j d\'A{ry,a)Q'^\r\a) (71) 
- j d^r'd\"A{ry ,Q.)FP^\r' y')p{r" ,a) 

with the ph-propagator 

Airy^O) = £ /^':"^^„ K.(r)K,(r-OK(r-)(^v2(rO (72) 

Vi,V2 Vl V2 

and the distribution of the transition strength 

S{a,A) = --Im f d\Q^^^{r,aYp{r,Q + iA) (73) 
n J 

These formulas refer to even-even nuclei because only in that case can one define the Landau 
quasi-particles as the single-particle excitations in the neighboring odd mass nuclei. This 
holds for spherical as well as for deformed nuclei lITSll . 
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We also have shown that the theory allows one to calculate moments and transitions in 
odd mass (A ± 1) nuclei, where A refers in general to a closed-shell nucleus. We then obtain 
for example, 

< A± l,a|e||3,A± 1 > - < A,0|e|0,A > 6„p = £ Q^%^P^f,, (74) 

V1V2 

with 

PviV9 ~ "ViaOy^p + T ~ 7^ ^ ^ViV3,V2V4Pv3V4 V'-'J 

c-vi tV2 '"^ap V3V4 
The derivation is given in |2|. Eqs. (l69l) and (1741) are the basis of the TFFS. 

3.5. Application to the giant electric dipole resonances in nuclei 

As examples for the TFFS and -later- for the ETFFS, we present numerical results for 
the giant electric dipole resonance in the doubly magic nuclei '**^Ca,^^^Sn and ^'^^Pb. The 
calculations presented here are performed in the r -space and not in the configuration space 
of a shell model. In the r -space one is able to include all single particle states which can 
contribute to a given multipolarity if the TFFS is applied, for example, to collective states 
in nuclei. The single particle continuum is included correctly using the method suggested 
by Shlomo and Bertsch lITlll . We wish to stress that the TFFS equation and the continuum 
RPA equation are formally identical. The TFFS equations are derived in a rigorous manner, 
therefore they are exact until one choses certain approximations for the quasi-particles, the 
interaction and the effective operators. As we have exact relations we are able to apply 
conservation laws, which connect the interaction parameters with gross properties of nuclei 
(symmetry energy, compression modulus and effective mass) and which give information 
about the renormalized operators (Ward identities). In addition, the Green function method is 
the most appropriate way to extend the theory in a consistent manner, as we will show in the 
next section. In the present calculations the experimental single particle energies are used as 
far as possible, and for the rest we take the solutions of a Woods-Saxon shell model potential, 
the parameters of which are adjusted to reproduce the experimentally known energies, as 
discussed in lUTI . 

For the calculation of the electric states we used the conventional Landau-Migdal 
parametrization for the ph-interaction: 

F''''(l,2)=Co5(ri-r2)[/o(p)+/^(p)x-x' + go(p)o-o' + g[)(p)o-o'x-x'] (76) 
with 

/(p)=/.x + (/,„-/«)poW (77) 
and with the parameters used in most of the previous applications. 

/„, = -0.002, U = -1-4, = 3,0, = 0.76, 

g = 0.05, g' = 0.96, Co = 300 MeVfm^ (78) 
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For the nuclear density Po('') in the interpolation formula we chose the theoretical ground 
state density distribution of the corresponding nucleus, 

PoW= E i(2jV + l)i?fW (79) 

which is more consistent than the previously used Woods-Saxon distribution. For that reason 
one had to readjust fgj^ and /^^ Here Ri{r) are the single-particle radial wave functions of the 
single-particle model used. 

For magnetic properties it has been shown lfT3l l6ll that one has also to consider the effects 
of the one-pion exchange. For that reason on must replace in the spin-isospin channel the 
conventional Landau-Migdal parametrization by the Julich-Stony Brook ansatz which has the 
following form in the momentum space: 

G'{q)= J -0^[V^{k)+Vp{k)]niq-k)+8G'oa-a'x-x' (80) 

with 

277-2 

Sl{q) = (27l)35(g) - — 5(|q| -q,) (81) 
T 

Here V^i and Vp denote the one-pion and one-rho exchange, t^(q) a parametrization of 
the nuclear correlation function [14J with qc = 3.93, which is the inverse of the Compton 
wavelength of the CO meson, and 5Gq a small (universal) correction parameter that must be 
adjusted to the data. The p exchange has been added because it introduces a natural cut-off 
for the tensor part of the one-pion exchange. A detailed discussion can be found in lfT3l l6l. 

The renormalized ph-interaction and therefore also the Landau-Migdal parameters are 
defined microscopically and can be traced back to the bare nucleon-nucleon interaction. 
In Ref.[6| these parameters have been calculated from a G-matrix based on a meson 
exchange model for the bare nucleon nucleon interaction. This is the lowest order for 
the renormalized ph-interaction nevertheless the theoretical results were in fair agreement 
with the phenomenological parameters. The Skyrme ansatz which is mostly used in self 
consistent nuclear structure calculations has the same form as the original Landau-Migdal 
parametrization. For that reason, the parameters of the Skyrme ansatz can directly be 
connected with the Landau-Migdal parameters These authors have shown that several 
parameterizations of the Skyrme forces led to instabilities against spin collapse. These 
relations can also be used to design new parameterizations for self consistent nuclear 
structure calculations [|5j|. For magnetic properties were the one-pion exchange has to be 
considered explicitly self consistent calculations based on chiral Lagrangian seems to be more 
appropriate. 

In the upper part of Fig.© the q-dependence of G' is plotted. One realizes that at 
small momentum transfers the central part of the spin-isospin interaction is strongly repulsive 
and for larger momentum transfers it is small. This is the reason why, for example, the 
Gamow-Teller resonances, where small momenta are involved, are strongly shifted to higher 
energies whereas high-spin states are not lfT3l l6ll. For the central part n -meson and p -meson 
contributions have the same sign and for the tensor part the opposite sign. The p -meson 
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Figure 2. Momentum dependence of the spin-isospin parameter G{q)' in units of 
[302MeVfm^] . The full line is the complete model ( I80> . the dashed line in the lower part 
is the correlated one-pion exchange only. 




exchange therefore acts as a natural cut-off for the strong tensor component of the one-pion 
exchange. 

We have already mentioned that the Migdal parameters have been calculated from a G- 
matrix starting from the bare nucleon nucleon interaction f6|. In this calculation the authors 
not only determined the zero-range parameters shows in eq. dTSb but also the momentum 
dependence of the spin-isospin interaction given in Fig.©. 

In the TFFS equations effective operators Q^^^ appear. In the case of electric multipole 
operators one can show using Ward identities that these operators are not renormalized due 
to the local gauge invariance, as shown in Appendix A. Magnetic operators are renormalized 
for two reasons: (I) there exists no conservation law for the spin contribution and (II) there 
are mesonic contributions. For that reason one has to parametrize the magnetic operators |l2l, 
the parameters of which are universal for all medium and heavy mass nuclei. In Fig.© the 
results of TFFS calculations for the electric giant dipole resonances in ^*^Ca, ^^^Sn and ^ospj-, 
are shown. The single particle continuum is included in an exact way. In order to simulate the 
experimental resolution and more complex configurations, a smearing parameter A =250 keV 
is used. The (isovector) electric dipole operator is given as: 

eN ^ eZ ^ 

Q\M = — Y, ^'^iM(a) - ^ E nYxu{^^) (82) 
1=1 i—\ 

In ^'^^Pb, where the electric giant dipole resonance is known ifTSll . the peak energies agrees 
with the experimental centroid-energies (up to 1 MeV) and the resonance exhaust a large 
fraction of the energy weighted sum rule. However, the theoretical width is much too small 
and the maximum of the cross section is much too high because the spreading width, which is 
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Figure 3. El photo absorption cross section for three double magic nuclei. The thin lines 
denote the results of the continuum RPA and the thick lines the results of the extended theory. 

due to the fragmentation of the single particle strength, is not considered. The same is true for 
the other two cases if one compares the theoretical results with neighboring nuclei data exist. 
In the next section we will discuss an extension of the TFFS, where the low lying collective 
phonons are coupled to the single particle states, which gives rise to a fragmentation of the 
single particle strength. In this way we include, as we will show, a major part of the spreading 
width. 

4. Extended theory of finite Fermi systems (ETFFS) 

As mentioned earlier, the original TFFS allows one to calculate only the centroid energies 
and total transition strength of giant resonances because the approach is restricted to Iplh 
configurations. In order to describe nuclear structure properties in more detail one has to 
include 2p2h, or even more complex configurations. A theoretical approach that takes into 
account the complete 2p2h configuration space and a realistic 2p2h interaction is numerically 
intractable if one uses a realistically large configuration space. For this reason the main 
approximation in the ETFFS concerns the selection of the most important 2p2h configurations. 
One knows from experiments, e.g. from the neighboring odd mass nuclei of ^^^Pb, that 
the coupling of low-lying collective phonons to the single particle states gives rise to a 
fragmentation of the single-particle strength, which is seen in even-even nuclei as a spreading 
width in the giant resonances. In the past 15 years an extension of the TFFS that includes, in 



Migdal Theory 



17 



a consistent microscopic way, the most collective low-lying phonons has been performed by 
Tselyaev lfT6ll and Kamerdzhiev, Speth and Tertychny ifTll (and references therein). In this 
way one considers a special class of 2p2h configurations because the phonons are calculated 
within the conventional (Iplh) TFFS. The phonons give rise to a modification of the particle 
and hole propagators, the ph-interaction and the ground state correlations. The extension that 
includes uncorrected 2p2h configurations is discussed in ref. [ 19 1 and L20il . 

4.1. Some basic relations 

In Landau's theory the self energy E is irreducible in the one-particle and one-hole channels, 
respectively, and the kernel K in the integral equation for the response function (see Fig.®) 
is irreducible in the particle-hole channel. One introduces now a hierarchy of energy 
dependencies: As one can see from eq. (l54b the particle-hole propagator introduces a strong 
dependence in the energy transfer Q.. Compared to this singular behavior one neglects in 
Landau's theory the energy dependence in the ph-interaction. One neglects consistently the 
energy dependence of the self-energy E in the Dyson equation (eq. (l34b and considers the 
quasi-particle and quasi-hole poles only. This approach is of leading order in the energy 
transfer Q.. In the extended theory, where phonons are introduced, one considers the next-to- 
leading order in the energy transfer; i.e., the self energy and the ph-interaction become energy 



dependent. We the write, explicitly, 

Ev(e)=Ev-FEf (£) (83) 

^lV3,V2V4(^) = ^iV3,V2V4 +^ViV3,V2V4(^) (84) 

> — S '''~>^ > > — S XT * 

« < <— <— — € — ^ « 

(a) (b) (c) (d) 



Figure 4. Graphs corresponding to the propagator of the (a) RPA and (b-d) the extended 
theory. The wavy Unes and the thin Unes denote the phonons and the single particle 
propagators, respectively. 

First we discuss the so-called g-^ approximation lfm i2ll. which is graphically shown in Fig.©. 
The upper part of graph (b) gives a correction to the Dyson equation for quasi-particles, which 
now has the form 

^vf (e) = 1 (85) 

Here gy'^'ic) is the one-quasi-particle Green function in the extended theory and the vertex 
7^!'', which couples the quasi-particle state v to the core excited configuration ij0i, is given by 

^^' = lCp5^aP (86) 



IV 1 



V2;i 
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where %' is the RPA wave function of the phonon considered. The corresponding energy- 
dependent correction for the ph-interaction has the form 



This model has been applied, for example, to Ml -resonances in closed shell nuclei ll2T1l . This 
approximation has, in general, problems with the so-called second order poles, which give 
rise to a distortion of the strength function near these poles. For this reason the model has 
been extended in several steps: (I) summation of the terms ll22t and (II) partial summation 
of diagrams |16|, which is call chronological decoupling of diagrams. In the latter case all 
lplh0 phonon contributions are consistently included and all more complex configurations 
are excluded -as long as one neglects ground- state correlations. The actual formulas, however, 
include the ground state correlations completely. In the present approach, like in lfT6l 
two types of ground state correlations are included: the conventional RPA ground state 
correlations, which affect the location and the magnitude of the residua of the ph propagator 
only, and the new type of ground state correlations, which is caused by the phonons. These 
correlations are qualitatively different from the conventional RPA correlations because they 
create new poles in the propagator, which then cause transitions between the Iplh (g) phonon 
in the ground state and the excited states. They give rise to a qualitative change of the strength 
distribution and a change in the sum rules for the moments of the strength function. As in the 
TFFS, the energy dependence of the generalized propagator A'^'^'^ is much stronger than that 
of the generalized ph interaction, as the next-to-leading order energy dependence is removed 
from the interaction and explicitly taken into account in the generalized propagator. Therefore 
one considers explicitly only the energy dependence of the propagator and neglects the energy 
dependence of the ph interaction. The interaction is parametrized as before; the corresponding 
parameters may differ from the previous ones. The final equation for the change of the quasi- 
particle density in an external field has the identical form to that in the TFFS ("eq. (l68t 



The analytical form of the generalized propagator can be found in Ref. lfTTl . 
4.2. Numerical applications 

As in the case of the TFFS the numerical input to the ETFFS are two sets of phenomenological 
parameters that describe (I) the Woods-Saxon single particle potential and (II) the effective ph 
interaction. If one includes phonons, the problem of double counting appears in the two input 
data sets. One therefore has to subtract explicitly the phonon contribution from the mean field 
solutions as well as the ph interaction. In order to obtain the new refined single particle basis 





(87) 



£ - £' + {Q-i - z6) 




(88) 
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^A, f solve the equation 



. 2 




(89) 



where the second term on the right side is graphically given by the upper part of graph 
(b) in Fig.®. As the wave functions are only little changed, we restrict the refinement 
to the single particle energies. In Fig.© we also show the ETFFS-results for "^^Ca, ^^^Sn 
and ^^^Pb. The effect of the phonon coupling is very clearly demonstrated in all three 
cases. The phonon coupling reduces the maximum of the strength distribution and shifts 
part of the strength to higher energies. This effect increases the width of the resonances 
and creates the long tails, which are seen in the data. For ^^^Pb the theoretical mean 
energy has the value Ef^ = 12,84MeV compared with experimental value Egxp = l3AMeV 
and the theoretical width Tth = 2>,6MeV compared to Texp = 4, IMeV. For The energy 
weighted sum rule (EWSR) is only little changed in the range between 5-30 MeV. For ^^Ca 
holds: Eth = \9,l\MeV compared to E^xp = \9,6MeV and Tjh = 5,36MeV compared to 
Texp = 7, IMeV . The phonon coupling improves the CRPA results by a factor of two. 

4.3. A realistic example for the application in data analysis 

The renewed interest in nuclear structure physics comes from newly built and planned 
experimental facilities that will give new information on nuclei far from stability. One of 
the main reasons for this is its application to astrophysics. Here, for example, one would like 
to know the compression modulus and the symmetry energy of neutron rich nuclei. Such 
experiments will be performed in so-called inverse kinematics; i.e. the radioactive isotopes 
will be accelerated and will be scattered from proton, deuteron or alpha particle targets. The 
analysis of such experiments will be very complicated as the breathing mode and giant electric 
dipole resonance are expected to be broad and because various multipole resonances overlap 
each other. For that reason on needs realistic strength distributions and energy dependent 
microscopic transition densities as input into the scattering codes in order to be able to extract 
reliable physical results. In the following we demonstrate this in the analysis of (a, a') cross 
sections in ^^Ni, where the isoscalar breathing mode has been investigated. A general problem 
in nuclei with A < 90 is that the isoscalar monopole resonance is very broad and no longer 
concentrated in one single peak. In deformed nuclei (such as one expects for nuclei far from 
stability) the giant resonances are split due the deformation, which gives rise to additional 
complications ll23l . In lighter nuclei the role of the surface becomes more important than 
in the heavier ones and, in addition, (a, a') experiments are most sensitive to the nuclear 
surface. For these reasons the correct treatment of the surface in those nuclei is crucial. All 
this has to be considered if one develops models for giant resonances in light and medium 
mass nuclei for which the conventional RPA approach is extended to include surface modes 
and the single-particle continuum. Compared to this extension of the configuration space, the 
way that one determines the mean field and the residual interactions is of less importance. 
When the present example was first published the experimental situation in ^^Ni was rather 
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Figure 5. Cross sections of ^^Ni(a,a') at 1.08°. At this angle the cross section is dominated 
by the isoscalar monopole resonance. 

unclear. In the original analysis of the (a, a') data j24j only 32% of the EWSR was observed. 
The authors used standard data analysis techniques with the same phenomenological transition 
densities for the whole energy range. For comparison the same type of experiment in 
Ca and a similar analysis by the same authors showed (92 ± 15) % of the EWSR ESl . 
With an improved -but still conventional- analysis the authors of Ref.|26| obtained about 50 
percent. Such differences might have serious consequences for nuclear matter compressibility 
and its application to astrophysics. In what follows we discuss an analysis of the (a, a') 
experiments in ^^Ni, where the isoscalar giant resonance has been analyzed using the results of 
ETFFS calculations. Within this model the authors calculated the distribution of the transition 
strengths of the isoscalar E0-E4 resonances and the corresponding energy dependent transition 
densities. Using these transition densities one obtains, in the standard way, (a, a') scattering 
cross sections that were compared with the data in the energy range between 5-35 MeV 
ll27l . One of the problems that occurs in such analysis is the subtraction of the instrumental 
background. Usually the experimentalists draw a straight line between the beginning and end 
of the resonance peak, where the strength distributions become flat. It has been demonstrated 
in Ref.[28| that if one adds up the contributions of the various weak multipole contributions 
one also obtains a flat back ground, which may erroneously be interpreted as part of the 
instrumental background. This background includes also a fraction of the resonance strength 
that one wishes to extract. In Fig.® the differential (a, a') cross section for the scattering 
angle 0=1 .08° is shown, where the EO contribution dominates the cross section. The full 
theoretical result (full line with dotes)[28| is compared with the experimental cross section 
ll27l . These data where obtained by subtracting an instrumental back ground denoted by 




Figure 6. Cross sections of ^^Ni(a,a') at 4.08°. At this angle the cross section is dominated 
by the isoscalar quadrupole resonance 



"Backgr(2000)" from the original experimental spectrum. The theoretical cross section is the 
sum of five different multipole resonances, of which the EO and E2 contributions are shown 
explicitly. The straight line in the lower part of the figure, denoted by " Backgr.(1996)" 
corresponds to the analysis of Ref.[24|, where the energy interval was smaller (12-25 MeV). 
In this analysis only 32% of the EWSR EO was found, whereas the present theory predicts 
71.4% of the EO EWSR in this regime. If one considers the larger 12-35 MeV interval, 
the theoretical strength exhaust 89.6%. We also show data taken at = 4.08°, where the 
EO resonance has a minimum and the E2 its maximum. The same analysis as before was 
performed by Ref . 11281 and the result is shown in Fig.®. Here the E2 strength dominates the 
cross section whereas the breathing mode contributes only a tiny fraction. The other isoscalar 
modes again give an important contribution to the cross section. As both experimental cross 
sections are well-reproduced within the ETFFS analysis one may conclude that the major 
part of the EO and E2 EWSR strength is in the energy range considered. One may also 
conclude that realistic microscopic nuclear structure models are necessary for the analysis of 
experiments in the giant resonance region. The models have to go beyond Iplh RPA and must 
include the coupling to low-lying phonons in order to reproduce the width of giant resonances 
in a realistic manner. 
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5. Conclusion 



We have reviewed Migdal's Theory of Finite Fermi Systems (TFFS) and the Extended Theory 
of Finite Fermi Systems (ETFFS) as well as their application to nuclear structure. Both 
approaches are derived within the Green function formalism, which enables one to obtain 
the final equations without approximations. For that reason one can apply conservations 
laws, which restrict the number of parameters in the actual numerical application. The 
TFFS and the ETFFS are frameworks that require single-particle energies, single-particle 
wave functions and the ph interactions as input. The single particle basis may be obtained 
from phenomenological single-particle models, as discuss here, or from self consistent mean 
field calculations. The original Landau Migdal ph-interaction includes delta functions and 
derivatives of delta functions in r -space, so that Skyrme forces can be directly connected 
with the Landau Migdal interaction [4J and L5|. This fact can be used as a test for the various 
versions of the Skyrme parameters but also for designing new set of parameters which are 
especially appropriate for self consistent nuclear structure calculations. The original Landau 
Migdal theory is formally identical with the Ip-lh RPA. Therefore if one calculates in this 
approach giant multipole resonances one can only reproduce the mean energy and the total 
strength of those states. For a quantitative description an extension of the theory is necessary. 
One has to include in a consistent way the effects of the phonons. This has been done in 
the Extended Theory of Finite Fermi Systems which we also reviewed here. We have pointed 
out that if one uses experimental single-particle energies and/or a correspondingly adjusted 
single particle model in the ETFFS, one has to correct for the phonon contributions in the 
single particle spectrum as the phonons (or, in general, 2p2h configurations) give rise to 
a compression of the single particle spectrum lfT9ll . Therefore a self consistent mean field 
approach provides the correct input data because in this case the phonon contributions are 
not included. For the same reason is the conventional RPA not the appropriate approximation 
in self consistent calculations because the uncorrected ph spectrum is in general too wide 
as a consequence of the phonon contribution being missing. In the case of very collective 
resonances, which are strongly shifted in energy, this deficiency can be removed by an 
appropriate ph interaction. W e are aware of the fact that, e.g. for a specific Skyrme 
parametrization the residual interaction is defined by the Skyrme interaction itself and one 
has, in principle, no freedom. But we are also aware of the fact that there exist dozens of 
different parameter sets that give similar ground state properties but quite different nuclear 
structure results Q. Moreover, the size of the configuration space can also be chosen in an 
appropriate way. Thus the arbitrariness in self-consistent calculations lies in the choice of the 
effective Lagrangian, effective Hamiltonian or energy-density functional and the size of the 
configuration space. For less collective states, such as the pygmy resonances, the energies of 
which are only little shifted by the residual interaction, self consistent calculations in general 
have difficulties reproducing simultaneously these experimental data and the giant electric 
dipole resonance. Therefore some authors, e.g.|29|, compress the too large level spacing by 
multiplying it adhoc with the corresponding effective mass in order to fit the data. This may 
be justified by the energy dependence of the effective mass but it certainly destroys the self 
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consistent treatment. We have also mentioned that for magnetic properties one has to consider 
in the spin-isospin channel the one pion exchange explicitly. For that reason self consistent 
calculations, starting from chiral Lagrangian may be the more consistent approach compared 
to the Skyrme Hamiltonian. 
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6. Appendix A: Ward identities 

As the final equations of the Landau -Migdal theory are derived without approximations, one 
can apply conservations laws in order to restrict the number of free parameters. Here we show, 
using local gauge invariance, that scalar vertices are not renormalized, i.e. one has to use the 
bare charges (for protons and neutrons) for the electric multipole operators. We start with the 
one-particle Green function: 

g{x,y) = -/ < AO|r (x) ^+ (y) } |A0 > (90) 

where x and y denote four- vectors. The gauge transformation of the ^-operators where q is 
the charge 

^''{x) = e''?^W^(x) ^'+(jc) = e-'^^W^+(x) (91) 
leads to the following change in the Hamiltonian 

m{x) = -nx)^ (92) 

Here f^^ = (j^p', j^p) is the four-current of neutrons and protons. This field does not give rise 
to any physical changes. If one expands the transformed operators in the one-particle Green 
function up to first order, one obtains 

^g{x,y) = i{qf{x)g{x,y) - g{x,y)qf{y)) (93) 

The inverse of this equation has the form 

5^"' = -g-'hg-' = -i{g'\x,y)qf{y)-qf{x)g-\x,y)) (94) 

On the other hand one can calculate the change of the Green function due to the additional 
term in the Hamiltonian eq. (l92b with help of the linear response function 

h{x,y) = j dzL{x,z;y,z)bH{z) (95) 
The corresponding inverse change of the Green function is given as: 

?>g-Hx,y) = J dzx{x,z;y,z)f{z)^ (96) 

With the field: 

f{z) = (97) 
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one obtains from eq. (l94l) and eq. (l96l) 



-i{g \x,y;(D + ^)qe'''y-qe'''''g \x,y-(D-^)) 



(98) 



iD. I d\x{x,z;y,z;G),Q.)p{z)e'''^ -i I d\x{x,z;y,z;(£),D.) {j{z)k)e 



Jkz 



This equation has in the single particle basis (py{x) the form: 

£ I (vi I ((0 + ^) IV3) (vake''^^ IV2) - (vi | qe''''' IV3) (V3 1 g"^ (w - y ) IV2) | 



(99) 



with 



and 



= ^2(vi|x(co,^2 
(Vi|T(co,i^; qe - 



qe 



ikr 



|v2) + (vi|?((0,^2; (j/:)e'^" )|V2) 



ikr 



)|V2) = £XviV3V2V4((0,^^) 
V3V4 



qe 



i{kr) 



J V4V3 



i{kr 



V4V3 



(100) 



(101) 



(vi|t(co,^2; {jk)e''''' ) |V2) = £ XviV3V2V4(o),^^) {jk)e'' 

V3V4 

We have replaced in eq.(l99t the charge distribution p by the point charge q. With the 
singular part of the one particle Green function eq. (l54l) one obtains a relation between the 
bare operators and the vertices x 



(0+§-evi (0-§-£v2 



(vi|<?e''"'|V2) 



Z\2 



(102) 
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We multiply this equation with ^/Zy^ and replace the vertices x by the renormalized vertices 
X defined in eq. (l60b 



{Q. - £vi - £v2) (Vl I {qe''"' |V2) = XviV2 



qe 



ikr 



+ X 



V1V2 



jke 



ikr 



(103) 



Here the approximation . ^ . was used. If one puts the renormalized vertex function 



Zv 



(eq AlOAl ) explicitly into eq. (fT03 

'^ViV3,V2V4 (W, f2) = X, 



CO 

'ViV3,V2V4 



(a),a) + 



(104) 



£ '^viv5,V2V6 ( ®' 
V5V6 V 



^V5 ^V6 
£v5 — £v6 ~ ^ 



^V5V3,V6V4 



£V5 + 



and considers the limit >> (£v[ — £v2)- One obtains the final result: 

EX,V2V4(^^'^''^)V3V4 = (^^'^'''^ 



;V3V4 ■ 



(105) 



V3V4 



Scalar (electric multipole) operators are not renormalized. 
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